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ABSTRACT 

In gamma-ray burst research it is quite common to fit the afterglow light curves with a broken 
power law to interpret the data. We apply this method to a computer simulated population of 
afterglows and find systematic differences between the known model parameters of the population 
and the ones derived from the power law fits. In general, the slope of the electron energy distribution 
is overestimated from the pre-break light curve slope while being underestimated from the post-break 
slope. We also find that the jet opening angle derived from the fits is overestimated in narrow jets 
and underestimated in wider ones. Results from fitting afterglow light curves with broken power laws 
must therefore be interpreted with caution since the uncertainties in the derived parameters might 
be larger than estimated from the fit. This may have implications for Hubble diagrams constructed 
using gamma-ray burst data. 

Subject headings: gamma rays: bursts — gamma rays: theory — methods: data analysis 



1. INTRODUCTION 

Many gamma-ray burst (GRB) afterglows can be in- 
terpr eted within the standard fireball model (see iPirai] 
120051 for a recent review), where a jet structure is implied 
from the steepening of the light curve ( Rhoads 1999). By 
fitting the light curves with a (sharp or smoothly joined) 
broken power law (see e.g. [Beuermann et al. 1999), it is 
possible to determine the time of the steepening (the so 
called jet break time) and the pre-break and post-break 
slopes of the light curves. This data can then be used to 
obtain information about the parameters of the underly- 
ing model, e.g. the jet opening angle which can then be 
used to est imate the energy of the burst (e.g. [F rail et al.' 
HlQI; Ghirla,nda et a,1.l l200l IFriedman B kK)m 2005^ . 

Although it has been shown that the afterglow light 
curves in the standard model can b e approximated with 
a broken power law (e.g. lSari et 311 11998: Rhoads 1999), 
the accuracy of the model parameters derived in that 
way has not been tested. In this letter we show that the 
results from fitting a broken power law to afterglow light 
curves must be interpreted with caution as the parame- 
ters derived from the fit can be systematically different 
from the actual model parameters. We do this by creat- 
ing a population of simulated afterglo w light curves using 
our v ersion of the standard model l|.T6hannesson et aLl 
l2005j) and then fit them with a broken power law. In 
section 2 we shortly describe the model and procedure 
used, present the results and then conclude in section 3 
with a discussion of our findings. 

2. THE MODEL AND RESULTS 

We used the instantaneous energ y injection version 
of the standard model presented in iJohannesson et al.l 
l|2005D . This is an extension of thelRhoadf^ 119991) model, 
and includes a more detailed calculation of the syn- 
chrotron emission and the effects of the equal arrival 
time surface (EATS). We created a population of 20,000 
afterglow light curves where the model parameter val- 
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ues were selected at random from a normal distribu- 
tion over a narrow parameter range given in table ^ 
We assume a constant density interstellar medium for 
all events. The standard deviation of the distribution 
for each parameter was fixed at 1/4 of its range and 
the distribution was clipped to fit within each param- 
eter range. The range of parameters was chosen after 
fitting several afterglows with our model o ver a wide 
range of frequencies (see examples i n Johanne sson et all 
l2005 Hde IJgarte Postigo et a l. 20051. Although the frac- 
tion of energy contained in electrons, eg, and magnetic 
field, tB, do not directly enter our formalism, they are 
included as they can affect the results vi a the chara cteris- 
tic synchrotron frequencies, Vm and Vc <)Sari "eTaLiri ggsi. 
The viewing a ngle of th e observer, 9^, can smooth the 
jet break ( Granot et al.i .2001'1 and was also included as 
a parameter. The light curves were evenly sampled in 
the logarithm of observer time and consisted of about 35 
points each (see fig. for a typical example). To make 
the synthetic data more realistic, we added normally dis- 
tributed random fluctuations with a standard deviation 
of 3%. We also assumed a fixed error of 5% for each data 
point. We find that our results are not sensitive to the 
values of these error parameters and lowering the error 
estimate only reduces the number of successful fits with- 
out reducing the scatter in our results. The results were 
also tested for the effects of the EATS by turning it off in 
a test sample. We found that it had no significant effect. 

According to standard fireball theory, the simplest af- 
terglow light curves can be described with the widely 
adopted analytical approximation ()Rhoads,.1999iL 
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for a given observing frequency, ly, in the range Vm < 
v < Vc- Here t is the time from the onset of the burst 
in the observer's frame, tj is the jet break time in the 
observer's frame and p is the electron energy distribution 
index. Since the jet break is in general not sharp, eq. |^ 
is often replaced with a smoothly joined broken power 
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TABLE 1 

The parameters of the afterglow population. 
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The parameters in the afterglow population arc normally 
distributed over the ranges shown here. The parameters are all 
distributed logarithmically, except p, 8v/6o and z. See text for 
more detail. 



law ijBeuermann et al]ll999|) . 
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Here, Fj is the flux at tj, —ai and —a2 are respectively 
the pre-break and post-break light curve slopes and n is 
a numerical factor controlling the sharpness of the break. 

The light curve break occurs when the jet enters the 
sideways expa nding regime (called exponential regime in 
IRhoadslll999|) and the observer receives light from the 
entire jet surface. The break time can be approximated 
as that point in time when 9 « 1/r, where 9 is the open- 
ing angle of the jet and F is the Lorentz factor of the 
relativistically moving shoc k front. Using a nalytical ap- 
proximations from Johanncs son et al.l lj2005j) we find that 



tj « 1.21(l+z) 
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where Eq is the total energy injected into the jet, no is 
the constant interstellar medium particle density, ^o is 
the initial opening angle of the jet in radians and z is the 
redshift of the burst. This formulation differs slightly 
from the one given in Rhoads (1999), because we choose 
to use the energy injected into the jet rather than the 
isotropic equivalent energy in order to better isolate 9o 
in the equation. 

As is often done when fitting GRB afterglow s (e.g. 
Wiiers et al.l E^997; Beuermann et al. 1999; Zeh et all 
we fitted the synthetic R-band light curves with 
both a sharply broken power law as in equation |^ 
and a smoothly joined broken power law as in equa- 
tion ^ . These will hereafter be referred to as sharp fits 
and smooth fits, respec t ively. The Levenberg-Marquardt 
method of IPress et al.l l|199fit) was used to minimize the 
value of the fit in both cases. Although the starting 
point of the fitting procedure was chosen as the theoreti- 
cally correct values of tj , ai , a2 and Fj , it did not result 
in an acceptable fit in every case. Only those events 
where the per degree of freedom is less than 1 for the 
smooth fits and 2 for the sharp fits were selected for fur- 
ther study. A higher threshold was used for the sharp 
fits since the numerically generated light curves are very 
smooth and not well represented by a sharply broken 
power law (see figure About half of the fits fulfilled 
those requirements in both cases. Since GRB afterglow 
measurements normally extend from a few hours after 
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Fig. 1. — A typical light curve from our sample. The points 
are from our model calculations, the solid curve shows the smooth 
fit and the dotted curve the sharp fit. The light curve break time 
determined by the smooth fit is rather accurate in this case, around 
22 days, whereas the sharp fit break time is around 3 days. The 
light curve slopes agree well with the data for the smooth fit but the 
post-break slope in the sharp fit is too shallow as may be expected 
from the break time error. 



TABLE 2 

The parameters of the Gaussian fits to the error 
distributions. 



Smooth Sharp 
Mean Std. dev. Mean Std. dev. 



Api/p 0.097 
Ap2/p -0.18 
Aej/Bo -0.034 



0.042 
0.064 
0.091* 



0.13 0.042 
-0.29 0.068 
-0.19 0.28** 



The Gaussian fit is not particularly good here and the distribu- 
tion is actually wider. 

There is a narrow peak around on top of the Gaussian. 



the burst to about a month, we limited our data and tj 
to this time range in the fitting procedure. Using equa- 
tion 0, it can be shown that this range also limits the 
range of derivable opening angles to about 1.5°- 15°, de- 
pending on other burst parameters. With real data, this 
range of opening angles can be further reduced if there 
is a bright underlying host or a supernova component 
making the afterglow light difficult to observe. 

The results for tj from the fitting procedure can be 
used to find an estimate of the opening angle of the jet 
by inverting equation (j2Jl, 
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Here, the opening angle is denoted by 9j to distinguish 
it from the known initial opening angle of the burst, ^oi 
from our numerical model calculation, although theoret- 
ically they should be equal. The results for the derived 
slopes, ai and a2, can be used to find the value of p and 
equation Q gives two different values, pi = (4/3)ai -I- 1 
and P2 — ct2- Comparison of the values obtained this way 
from the fitting procedure to the parameters known from 
the numerical model is shown in figure El as the distribu- 
tion of relative differences between the derived param- 
eters and the known model parameters: Api/p, Ap2/p 
and A9j/9Q, where Api = Pi — p, Ap2 = P2 ~ P and 
A9j — 9j — 9q. Those distributions were each fitted with 
a Gaussian profile that is overlaid on the distributions in 



3 




-0.6 -0.4 -0.2 -0.0 0.2 0.4 0.6 
Ap..2/p 



2000 




Ae/e„ 

Fig. 2. — The distribution of relative differences in the value 
of parameters deduced from the fits. The upper panel shows the 
distribution in Api/p and Ap2/p- The smooth and sharp fits are 
represented, respectively, by plus-signs and diamonds for Api/p 
and triangles and boxes for Ap2/p- The lower panel shows the dis- 
tribution in A8j/8o for smooth (plus-signs) and sharp (diamonds) 
fits. The overlaid curves are Gaussian fits to the corresponding dis- 
tributions. The parameters of the Gaussians are given in table 121 
The count is higher in the lower panel because the binsize is larger. 



the figure. The parameters of the Gaussian profiles are 
presented in table [3 

The distributions and the corresponding Gaussian fits 
clearly show a systematic difference in the evaluation of p, 
where pi is overestimated and p2 underestimated. From 
the A6j 1 9o distribution, it is clear that the sharply bro- 
ken power law does not do a good job in determining 
the opening angle and there is both a significant under- 
estimate and a large standard deviation. This is also re- 
flected in the Ap2 /p distribution since an underestimate 
of the opening angle will make p2 smaller. It should how- 
ever be noted that there is a high narrow peak around 
zero in the A9j/9o distribution for the sharp fits which 
is not included in the Gaussian fit. The overall stan- 
dard deviations in the Gaussian fits are not very high 
but should be considered when using 9j to correct the 
isotropic energy, since the correction factor is approxi- 
mately proportional to 

We have checked for correlations between the parame- 
ters pi, p2 and 9j deduced from the fits and the known 
parameters of our numerical model, which range is shown 
in tabled We find no significant correlation within that 
parameter range, except between 9q and A9j/9o. Fig- 
ure O shows a scatter plot of A9j /^o as a function of 
for both smooth and sharp fits and also a binned average 
shown as a solid line. It is clear that narrow jet opening 
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Fig. 3. — The correlation between 9q and A9j/6o for smooth 
fits (top) and sharp fits (bottom). In these scatter plots, the solid 
curves show a binned average of the data. Also shown is the binned 
average calculated from tj^Y (dashed curve). 



angles are systematically overestimated and wider ones 
underestimated. This is more evident with the sharp 
than the smooth fits. When the jet opening angle is used 
to correct the isotropic energy for beaming, this will lead 
to a clustering of the derived energy values. 

To check the validity of the approximations used in 
deriving equations l|2Jl and I0J we determined the char- 
acteristic time tj^r when 9 equals 1 /F in our model calcu- 
lations. We then estimated the opening angle by setting 
tj equal to tj^r in equation Q). This should give correct 
results in each case if these approximations are valid. 
The resulting angle, 9j, is then compared to the known 
model opening angle, 9q, in a similar way as the for the 
power law fits. Figure 13 also shows the binned average 
of the correlation between A9j/9o and 9o for 9j calcu- 
lated from tj.r (dashed curve). It clearly shows that the 
approximations used in deriving equations © and 
break down for very narrow jets. This is because the 
j et break time can be ve ry close to the deceleration time 
ijPanaitescu et al.lll998j) where the approximations are 
not valid. These approximations also break down in very 
wide jets since these may not necessarily be assumed to 
be ultra relativistic at the time of the jet break. Fig- 
ure |3 also shows that the correlation is stronger when 9j 
is calculated from tj (solid curves), so tj is overestimated 
from the fits when the jet is narrow and underestimated 
in wide ones. One explanation is that the jet break time 
is not necessarily within the limited time range used in 
the fits, but it could also be due to j/,„ crossing the ob- 
serving frequency, especially in narrow jets. In the latter 
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case, those events can be identified from a small value of 
the smoothness factor, n, because the break around Vm 
is much smoother than the jet break. This crossing of Vm 
through the observing frequency is also the cause of the 
small wings in the distribution of Api/p and Ap2/p seen 
around the value -0.1 in figureEl By looking at the indi- 
vidual spectra, we can identify those events where Vm or 
Vc cross the observing frequency. We find that the latter 
case does not have a significant effect on our results. 

3. CONCLUSION 

We have shown that results from the standard pro- 
cedure of fitting afterglow light curves with a broken 
power law must be interpreted with caution. There can 
be systematic differences in the evaluation of the elec- 
tron energy distribution index from the slope of the fit- 
ted curves and a strong correlation between the relative 
difference in the opening angle estimated from the light 
curve break time and the initial opening angle of the 
jet. These findings are partly due to the approximations 
used in deriving equations (l)-(4) being used out of their 
validity limits. This applies to appro ximations in both 
the dynamical (|Bianco fc Ruf5nill2005 D and the radiation 
properties of the expanding shell. These differences and 
particularly the correlation is also a consequence of diffi- 
culties in accurately determining the jet break time from 
the afterglow light curves. 

The fitting procedure we used was completely auto- 
matic and often did not converge. In some cases it 
resulted in erroneous parameter values and we there- 
fore adopted a threshold on the reduced to eliminate 
those bad fits. To test the effect of this threshold, we 
removed it from the selection criteria and re-examined 
the data, still considering only those events where the 
light curve break time was within the time range of our 
data. This left us with over 90% of the original popula- 
tion. The most significant changes we find in the results, 
were slightly larger standard deviations in the relative 
difference distributions and more spread in the A9j/9Q 
- 9q correlation thereby weakening it. The spike around 
seen in the A6j/6o distribution for the sharp fits in 
figure El (diamonds) also becomes a dominating feature. 
This indicates that by using the as the strongest filter, 
we remove the bad fits from the ensemble but may also 
lose some useful fits. 

It is known theoretically that the magnitude of the 
jet break, Aa = 02 — ai, can be used to differentiate 
bet ween a wind like or a c onstant density environment 
fe.g-. iPanaitescu et alJll998 l). The systematic differences 
we find in the evaluation of p indicate that Aa is under- 



estimated in most of our events. Similar conclusion is 
obtained for a population of afterglow light curves in a 
wind medium, the main difference being that pi did not 
show a systematic deviation in that case. This renders 
the method of using Aa to distinguish between density 
profiles impractical. 

The correlation between the opening angle estimate 
and the known opening angle puts a strong limit 
on interpretati ons of the beaming corrected energy of 
the burst (e.g. Frail ct al. 2 001t iGhirlanda et all l2004t 
lFriedma~fc Blo om 2005). The clustering of the jet 
break time due to limited time span of the data together 
with a generous use of the approximations used in de- 
termining Oj results in a bias towards moderate opening 
angles, approximately between 2°- 10°. The rather large 
standard deviation in the A9j/9o distribution also makes 
the results unrehable. Using the beaming corrected en- 
ergy as a basis for cosmological studies therefore calls for 
a very careful determination of tj and 9j . 

The synthetic burst population studied in this paper 
is computer generated and the power law fit should in 
theory be perfect. The fact that the differences be- 
tween the parameters derived from the fits and the known 
model parameters are so significant, makes the accuracy 
of power law fits to real measurements of afterglow light 
curves a concern. In real bursts, effects such as den- 
sity fluctuations and energy injection, can change the 
shape of afterglow light curves as may for e xample be 
the case for GR B 021004 Ce.g. iLazzati et fd\ \2(m] and 
GRB 030329 fe.g. lSheth et a l.''2003'). These were densely 
sampled and were not well fitted with a broken power 
law due to bumps in the light curves. It is not hard 
to see that grainier measurements of the same events 
could have been fitted with a broken power law, leading 
to even larger uncertainties in the parameter estimates 
than discussed here. It should also be noted that other 
models are capable of explaining light curve breaks, the 
most popular being the structured jet model ijRossi et alJ 
1200^ . There the fight curve break depends on the ob- 
server's viewing angle rather than the jet opening angle. 
Hence the methodology adopted in this letter is not di- 
rectly applicable to that model. 
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